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• ^ ! 1 Introduction 

X: 

H , The fractional Brownian motion (fBm) with Hurst parameter H G (0, 1) is a 

. - - i zero mean Gaussian process = {B^,t > 0} with covariance fmiction 



RH{s,t) = l[t'" + s'^-\t~s\'"). (1.1) 



This process was introduced by Kolmogorov in [12] and later studied by Man- 
delbrot and Van Ness in [17]. Its self-similar and long-range dependence (if 
H > ^) properties make this process a useful driving noise in models arising in 
physics, telecommunication networks, finance and other fields. For a review of 
some applications of fBm we refer to [3] . 
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From (jl.ip we deduce that E( I — ) = |i— and, as a consequence, 
the trajectories of are almost surely locally a-H61der continuous for all 
a G {0,H). Since B^ is not a semimartingale if H ^ 1/2 (see [10]), we cannot 
use the classical Ito theory to construct a stochastic calculus with respect to the 
fBm. Over the last years some new techniques have been developed in order to 
define stochastic integrals with respect to fBm. Essentially two different types 
of integrals can be defined: 

i) The Skorokhod integral (or divergence integral) with respect to fBm is 
defined as the adjoint of the derivative operator in the framework of the 
Malliavin calculus. This approach was introduced by Decreusefond and 
Ustiinel [5 , and later developed by Carmona and Coutin [2], Duncan, Hu 
and Pasik-Duncan [7], Alos, Mazet and Nualart [T], Hu and 0ksendal [TT] . 
Cheridito and Nualart [31, among others. This stochastic integral satisfies 
the zero mean property and it can be expressed as the limit of Riemann 
sums defined using Wick products. 

ii) The pathwise Riemann-Stieltjes integral VgdB^ exists if {vt,t G [0, T]} 
is a stochastic process with Holder continuous paths of order a > 1 — H, 
as a consequence of the results of Young [51] . Zahle expressed this integral 
in terms of fractional derivative operators, using the fractional integration 
by parts formula (see [H]). We also refer to [5] for a pathwise approach 
to the stochastic calculus based on the regularization of the noise. 

The aim of this paper is to study the d-dimensional stochastic differential 
equation 

Xt=Xo+ [ b{s,X,)ds+ [ aw{s,X,)dWs+ [ aff(s,X,)dSf, (1.2) 
Jo Jo Jo 

where W is an r-dimensional standard Brownian motion and B^ is an m- 
dimensional fractional Brownian motion with H e (51 1)- We assume that 
the processes W and B^ are independent. In the above stochastic differential 
equation, the integral aw{s, Xs)dWs should be interpreted as an Ito stochas- 
tic integral and the integral J* anis, Xs)dBs as a pathwise Riemann-Stieltjes 
integral in the sense of Zahle [2B . Our main result is a general theorem about 
the existence and uniqueness of solutions for the stochastic differential equation 
()1.2|) under suitable conditions on the coefficients. 

Equations driven only by a fBm with Hurst parameter H G (5,1) can be 
solved by a pathwise approach using the p- variation norm (see [15j). the frac- 
tional calculus introduced by Zahle (see [35] and [H]), or Holder norms [3T] . 
Also using the tools of rough path analysis introduced by Lyons in [16j , Coutin 
and Qian proved in ^ the existence of strong solutions for stochastic differential 
equations driven by fBm with H > ■j and studied a Wong-Zakai approximation 
limit for these stochastic differential equations. 

Kubilius has studied stochastic differential equations driven by both fBm 
and standard Brownian motion (see [T3]), in the one dimensional case, with 
<^w J c// independent of the time and with no drift term (6 = 0). In this setting, 
the author proves an existence and uniqueness result provided that aw is a 
Lipschitz function and (Th G C^~^^ , with 6 > q{l — H), q > 2. With these 
assumptions, the solution belongs to the space of continuous functions with 
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g-bounded variation. Kubilius defines the stochastic integral with respect to 
fBm as an extended Riemann-Stieltjes pathwise integral and he uses p- variation 
estimates. 

Our approach is completely different from Kubilius [13] in the sense that 
we combine the pathwise approach with the Ito stochastic calculus in order to 
handle both types of integrals. Then, the uniqueness of a solution follows from 
estimates for both Ito and Riemann-Stieltjes integrals. However, the existence 
of a strong solution cannot be obtained by the classical fixed point argument 
because the estimates of the Holder norm of an integral with respect to 
produce some higher order terms. For this reason, we first prove the existence 
of weak solutions and later deduce the existence of strong solutions using the 
Yamada-Watanabe theorem. 

The paper is organized as follows. In Section 2 we state the problem and 
list our assumptions on the coefficients of Eq. ()1.2p . Section 3 provides some 
estimates for fractional and Ito integrals. In section 4, the pathwise uniqueness 
property of solutions of Eq. (jl.2p is proved. In section 5, we introduce the Euler 
sequence that approximates the solution of the stochastic differential equation 
and prove that it is a tight sequence. Then, the Skorokhod representation 
theorem is applied in order to prove the existence of a weak solution for the 
stochastic differential equation. Finally, we prove the existence of a unique 
strong solution by using the Yamada-Watanabe theorem. 

2 Preliminaries 

Fix a time interval [0,T] and a complete probability space (O, J-", P). Suppose 
that = {Bl^,t G [0, r]} is an m-dimensional fractional Brownian motion 
with Hurst parameter H G and W = {Wt,t G [0, T]} an r-dimensional 

standard Brownian motion, independent of B^ . Let Xq be a d-dimensional 
random variable independent of {B^ , W). For each t G [0, T] we denote by Tt 
the cr-field generated by the random variables {Xq, B^ , Ws, s G [0,t]} and the 
P-nuU sets. 

In addition to the natural filtration {Ttjt G [0,r]} we will consider bigger 
filtrations {Gt^t G [0,r]} such that: 

1- {Gt} is right-continuous and Go contains the P-nuU sets. 

2. Xq and B^ are ^/o-nieasurable, and is a tJt-Brownian motion. 

Notice that J-t C Gt, where is the u-field generated by the random vari- 
ables {Xo,B", Ws,s G [0,t]} and the P-null sets. 

Consider the stochastic differential equation (|1.2p , where Xq is a d-dimcnsional 
random variable independent of {B^ ,W) and the coefhcients are measurable 
functions 6* , ct^^ , cr}^^' : [0, T] x M'' ^ M, 1 < i < d, 1 < A: < r, 1 < j < to. We 
will make use of the following assumptions on the coefhcients. 

(H6) The function b(t, x) is continuous. Moreover, it is Lipschitz continuous in 
the variable x and has linear growth in the same variable, uniformly in t, 
that is, there exist constants Li and L2 such that 

\b{t,x)-b{t,y)\<Li\x-y\, 
\bit,x)\<L2{l + \x\), 
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for all x,y eR'^ and t e [0, T]. 

(HiTvi') The function aw{t,x) is continuous. Moreover, it is Lipschitz continuous 
in X and has linear growth in the same variable, uniformly in t, that is, 
there exist constants L3 and L4 such that 

\aw{t,x) - aw{t,y)\ < is k - ?/| , 
\<Jw{t,x)\ <L4(l + |x|), 

for aU x,y eR'^ and t e [0, T]. 

(HcTff ) The function anit, x) is continuous and continuously differentiable in the 
variable x. Moreover, there exist constants ^5, Lg and L7 such that 

\dxi<TH{t,x)\ < i5, 

\dj;^aH{t:X) - d.j,^aH{t,v)\ < Lq\x - yf , 
\aH{t,x) - (Th{s,x)\ + \dxi<JH{t,x) ~ dx,cTH{s,x)\ <Lj\t- sf , 

for all X, J/ e M'' and t G [0, T], and for some constants < 5,(3 < 1. 

Note that assumption (Her//) implies the linear growth property, i. e., there 
exists a constant L such that 

\aH{t.x)\<L{l + \x\) (2.1) 

for aU x,y(^R'^ and t £ [0, T]. 

Let us now introduce some function spaces that will be used in the analysis 
of solutions of the stochastic differential equation p.2p . Let < a < ^. For 
any measurable function / : [0, T] — > R"^ we introduce the following notation 

ll/WIL:=l/WI+ r '^f ~{l?' d^- (2-2) 
Jo (i — s) 

Denote by Wq'°° the space of measurable functions / : [0, T] — > R'^ such that 

ll/IL,oo:= sup ||/(t)IL<«). (2.3) 

te[o,T] 

For ^ e (0, 1] let C be space of /x-Holder continuous functions / : [0,r] R'', 
equipped with the norm 

11/11^:= II/II00+ Bup ttl/(!^<oo, (2.4) 

0<s<t<T [t~S) 

where := sup \f{t)\ . Given any e such that < e < a, we have the 

te[o,T] 

following inclusions: 

C"+' c M^o"'°° C C"-'. (2.5) 

In particular, both the fractional Brownian motion , with H > ^, and the 
standard Brownian motion W, have their trajectories in Wq'°°. We denote by 
the conditional expectation given Tq, that is, given Xq and . 
We now define the space of processes where we will search for solutions of 
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Definition 2.1 LetWg be the space of d- dimensional Qt-o-dapted stochastic pro- 
cesses X = {Xt,t G [0,r]} such that almost surely the trajectories of X belong 



to Wo'°" andJ^E^ \\X, 



ds < oo. 



A strong solution of the stochastic differential equation (|1.2p is a stochastic 
process X in the space Wjr, which satisfies (ll.2p a.s. The main result proved in 
this paper is the following theorem on the uniqueness and existence of strong 
solutions for (11.21). 



Theorem 2.2 Assume that the coefficients 6, crw and an satisfy the assump- 
tions (Hb), (Haw) and (Han)- If ot satisfies 1 — H < a < min{i,/3, |}, then 
there exists a unique strong solution X of Equation 

Remark Notice that in all our results we can replace the fractional Brownian 
motion by an arbitrary stochastic process with Holder continuous trajecto- 
ries of order 7 > ^. 

3 Integral estimates 

In this section we will first define the integral with respect to fBm as a gener- 
alized Stieltjes integral, following the work of Zahle [26]. We also present some 
basic estimates of this integral. 

Let / e L^{a,b) and a > 0. The left-sided and right-sided fractional 
Riemann-Liouville integrals of / of order a are defined for almost all x G (a, b) 

by 

1 



and 



i:+fi^)--=7^^ I {x~yr-'f{y)dy 
r(a) Ja 



1 



ib-fi^)--^ff^^ I iy-xr-'f{y)dy 



respectively, where T{a) := ^" ^dr is the Euler gamma function. Let 
/"+ {LP) (resp. (L^)) be the image of LP{a, b) by the operator /"+ (resp. ). 
If / e (LP) (resp. / e 4" (LP)) and < a < 1 then the Weyl derivatives of 
/ are given by 

^''^^(")^=f(T3^(,(^3^ + "X (,-,).+! ^^j W)W (3-1) 

and 

r>a ft 1 ( fix) , f'' fix) - fjy) ^ \ 1 

^^-^(")-=f(r3^1^(6^+"X 7^3^dyjl(.,,)(.), (3.2) 

respectively, and are defined for almost all x £ (a, b) (the convergence of the 
integrals at the singularity y ~ x holds pointwise for almost all x G (a, b) if 
p = 1 and moreover in L^'-sense if 1 < p < 00). 
We have that: 

. If a < i and <? = ^ then /^"+(Lf ) = (LP) C i«(a, b). 
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• If a > i then {LP) U /^"_ (LP) C C""p (a, b). 

The fractional integrals and derivatives are related by the inversion formulas 

and similar formulas also hold for and D^^ . We refer to [H] for a detailed 
account on the properties of fractional operators. 

Let f{a+) := lim/(a + e) and g{h~) := limg(6 — e) (we are assuming that 

£\0 e\0' 

these limits exist and are finite) and define 

:= (/(x)-/(a+))l(,,,)(a;), 
9b-{x) ■■= {g{x) - g{b~)) l{a,b){x)- 

We recall from the definition of generalized Stieltjes fractional integral with 
respect to irregular functions. 

Definition 3.1 (Generalized Stieltjes Integral) Suppose that f and g are 
functions such that f(a+),g{a+) and g{b—) exist, /„+ S I^^(LP) and gi,- G 
I^Z'^{LP) for some p,q> 1/q < 1,0 < a < 1. Then the integral of f 

with respect to g is defined by 

f fdg := (-1)" /' D:^f{x)Dll"g,- ix)dx + f{a+) (5(6-) - <?(a+)) . 

J a J a 

Remark 3.2 The above definition is simpler in the following cases. 

• If ap <\, under the assumptions of the preceding definition, we have that 
f S {L^) and we can write 



f fdg=[-ir f D:^f,Ax)Dll" 



9b-(x)dx. (3.3) 

• If f IE C^{a,b) and g G C{a,b) with A + > 1 then (see ^1) we can 
choose a such that 1 — < a < A, the generalized Stieltjes integral exists, 
it is given by i3.S\) and coincides with the Riemann- Stieltjes integral. 

The linear spaces (LP) are Banach spaces with respect to the norms 

ll/ll/»^(L.) ■■= II/IIl. + Pa./llL. - P^/IIl. > 

and the same is true for the spaces I^^{LP). If ap < 1 then the norms on 
/"+ [LP) and {LP) are equivalent and \i a < c < d <b, then 

fdg := / l(c,d)fdg. 

J a 

Now, fix the parameter a such that < a < ^, denote by the space 

of measurable functions g : [0, T] R™ such that 

II II f \9it)-9{s)\ , f \g{y)-g{s)\ ^ \ ^ 
ll5lli-a,oo,T — sup ( — + / — r^^jTT^dy ) < 00. 



0<s<t<T 



(t - s)i-" J, {y - s) 
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and denote by W^'^ihe space of measurable functions / : [0, T] M.'^ such that 



l/WI 



ds 



I ^ — r;— dyds < oo. 

Jo is-yr+' 



It is easy to prove that 



for all e>0. For g G Wl we have that 



Aa(5) 



1 



< 



r(l - a)o<s<t<T 
1 



sup |(Z?i-"gt_)(s)| 



\l-a,OD,T 



< OO. 



r(l -a)r(a) 

Moreover, if / G Wq''^ and g G then fdg exists for all t G [0, T] and 



/d.g 



<Aa(3)||/La 



(3.4) 



Now we will deduce useful estimates for the integrals involved in Equation 
()1.2p . Fix a G {1 — H, i). We will denote by C a generic constant which depends 
on the constants Li, I < i < 7, (3 and S in the assumptions, on T, a and the 
dimensions r,d,m. For any function / G Wq'"^ define 



Jo 

Proposition 3.3 If f e Wo'°° then (/) G W^'°° and for all t G [0, T] 



\Fhf)\\^<c(^f^ 



1/(^)1 



ds + 1 



(3.5) 



Proof. By Proposition 4.3 in [1^] and the growth assumption in (Hfe) we 
have that 



\F^{f)t<C f 
Jo 



\b{s,fis))\ 
(t - s)" 



ds 



<cr\im^ds 



< c 



it~s)" 
* 1/(^)1 



(t - s)" 



ds + 1 



Proposition 3.4 // /, h G W^o"^°° f/ien /or all t G [0, T] 



\F:^if)-F^ih)\\<C 



' ll/(^)"M^)L 

(t-s)" 



ds. 



(3.6) 
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Proof. By Proposition 4.3 in [19^ and the Lipschitz assumption in (H6), we 
have that 

||f,V)-^.>)IL<c/M^JM,, 

Jo [t - s) 

■ 

Given a function / £ Wq'°°, let us define 

G^a) := f aH{sJ{s))dBf. 
Jo 

Proposition 3.5 Suppose 1 — H < a < min(i, /3). Then for all t G [0, T] 

\\G^,"{f)\l<CK{B") f {{t-s)-^- +,-)(i + ||/(,)||jd.. (3.7) 



Proof. By Proposition 4.1 of 19j and the Holder continuity in time, given 
in assumption (Hcr/f ), we have 

Jo 

<CK{B") [\{t-s)-'" +s-")il + \\f{s)\Uds. 



Proposition 3.6 If 1 — H < a < min and f , h £ Wq'°° then for all 

te [o,r] 

\\Gl"{f)-Gl''{h)\\^<Ck^{B") (3.8) 

X f{{t-s)-^'^ +.-")(l + A/(.) + AM.))||/(s)-M«)Lds, 
Jo 

where we denote 

A,(„..rM#a.. (3.8, 

Jo (s - r) 

Proof. From Proposition 4.1. of T^, we have that 

Jo 

X hnis, f{s)) - crnis, h{s))\\^ ds 

<CKo.{B") f {{t-s)'^" +5-")(|aff(s,/(s))-aH(s,/i(s))| 
Jo 

\oh{s, f{s)) - (Tff(s, h{s)) - anir, f(r)) + auir, h{r))\ 



/ X , 1 dr ds. 

(s - r)«+i 



Now, using the assumptions in (Hcr^f) and Lemma 7.1 in [19j . we have that 

\a{ti,Xi) - (T{t2,X2) - <T{ti,Xz) + (7{t2,Xi)\ < 
< L5 \xi - X2 - X3 + X4: \ + L-! \xi - X'i \ |<2 - tif 
+ Lg \xi - X^\ {\xi - X2t + - Xi] 

As a consequence, 



\f{s)-h{s)-f{T) + h{r)\ 



(s - r)" 



+1 



dr ds 



Jo 

X (1 + (A/(,s))^ + {Ah{.s)f) ||/(.) - his)\\l ds. 
■ 

Finally, we will consider the Ito stochastic integral with respect to the r- 
dimensional standard Brownian motion W. The following lemma is an imme- 
diate consequence of Ito calculus. 

Lemma 3.7 Suppose that u = {u{t),t e [0,T]} is an r- dimensional Qt-adapted 
stochastic process such that /j, [u(sf] ds < 00. Then for all t G [0, T] a.e. 

(3.10) 









2" 








/ u{s)dWs 










Jo 


a 


Jo 



< C / (i-s)-^-"E^ ds 



Proof. Notice first that, by Fubini's theorem, the right-hand side of p.lOp 
is finite for all t G [0, T] a.e. Applying the Ito isometry property and the 
Cauchy-Schwarz inequality, we have that: 

21 



E" 



< C 







2" 












f u{s)dW, 




< CE^ 




/ u{s)dWs 


■•( 




Jo 


a 






Jo 





/; u{r)dWr 



[m(s)2] ds+ - 



T5- 



t E 



w 



J' uir)dWr 



Jo {t- s) 



(t-s) 



-ds 



a+1 



ds 



Therefore, by the Ito isometry and Fubini's theorem we obtain 



w 



E 



< C 



< c 



u{s)dWs 



u(s) J ds -f- yj 



E^ [u{sY\ ds + jj 



2 -'O 



(<-s)"^""ds 



E^ [u{rf] dr 



E"' [u{rf 



-dr 



<C (i-O^^-^EVV \^u{rf] dr. 
Jo 
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Assume that / = {f{t),t £ [0,T]} is a d-dimensional stochastic process in 
Wg. Define 

era) := f <j{sJ{s))AWs. 
Jo 

We have the following estimates for these integrals: 
Proposition 3.8 Let / £ Wg. Then for all t e [0, T] a.e. 



E 



w 



Gr(/) L <c (t-s) 



1 



II/WII' 



ds. 



(3.11) 



Proof. It follows from (j3.10p and the linear growth assumption in (Hcrvi"^ 



Proposition 3.9 Let f,h eWg. Then for all t e [0,T] a.e. 



E 



w 



Gnf)-G? 



w 



<C I (t-s)-^-"E'^ \f{s)~h{s)\' 

'0 



ds. 
(3.12) 



Proof. By estimate p.lOp and the Lipschitz assumption in (Raw), we 
obtain 



E 



w 



or if) -or mi 



<C (t-s)- 
Jo 

<C At-,s)-^-"E^ \\f{s)-h{s)\^ 



e^wis, /(s)) - crwis, h{s))\ 
ds. 



ds 



4 Pathwise uniqueness 



In this section we define the notion of weak solution for the stochastic differential 
equation (|1.2p and we discuss the pathwise uniqueness of a solution. 



Definition 4.1 A weak solution of the stochastic differential equation is 
a triple {X,B",W) , {n,T,P), {Gt^t e [0,r]}, where 

1. {Cl,J-,P) is a complete probability space, {Gt} is a right- continuous fil- 
tration such that Q{) contains the P- null sets. 

2. W is a Qt-r- dimensional Brownian motion. 

3. is a fractional Brownian motion of Hurst parameter H which is Qq- 
measurable. 

4-. The process X is Qt-adapted, has trajectories in W^'"^ almost surely, and 



Jo 



ds < cxD a.s. 
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5. (X,B^,W) satisfies Equation U.2\) a.s. 

Definition 4.2 We say that pathwise uniqueness holds for Equation if, 
whenever (^X, W, B^^ and (Y, W, B^^ are two weak solutions, defined on the 
same probability space (Jl, J- , P) with the same filtration {Gt} o-nd Xq — Yq a.s., 
then X ^Y. 

We will make use of the following technical lemma. 

Lemma 4.3 LetO < r/ < 1/2. /// is a continuous function such that \\f\\^ < N 
and a < rjS, then A(/) is bounded by a constant C depending on T , N , a, 5, 
and rj, where we use the notation introduced in \2.4^ and h3.9^) . 

Proof. Clearly 

A(/)M^ri^w-/':ii'd.-<iv'^, 

^ ' - ' Jo (s-r)"+i - 77(5 -a' 

which gives the result. ■ 

Let / G Wq'°°. By the estimates proved in Proposition 4.2 and Proposition 
4.4 of [H], the sample paths of the integral processes F^{f) and G'^^ [f) are 
continuously differentiable and Ty-Holder continuous of order 77 < 1 — a, respec- 
tively. Therefore, if X is a weak solution of ()1.2|) . then the trajectories of X are 
yy-Holder continuous for all 77 < 1/2. 

Theorem 4.4 (Pathwise uniqueness) Let 1 — H < a < min{/3, 



Then, the pathwise uniqueness property holds for Equation U.S^) . 

Proof. Let X and Y be two weak solutions of ()1.2|) defined on the same 
probability space, adapted to the same filtration and with the same initial value. 
Then the trajectories of X and Y are ry-Holder continuous, for all rj < 1/2. 
Choose rj such that a < j] < 1/2. Consider the sets ft^ C fl, defined by 

nN:={Luen: \\X\\^ < N and < TV } , 

with N E 'N. It is clear that f^Ar fi. From (|1.2p we have that the difference 
between the two solutions satisfies 



< 4 



\\Xt-YtCln, 
:^\\{Fl'{X)~F^{Y)) 1, 



-4E 



w 



(Grw-Gr(n)iil 



-4E 



w 



{G^"iX)-G^"{Y))U 



(4.1) 



We split the set ft into il^ and il\ilN in the second summand of (|4.ip and use 
the estimates p.6p . p.8p . (|3.12p in order to obtain 



w 



E 



< C 



\Xt~Y\ 



\Xf, — Ys 



^{s,t)E^ 
ip{s,t) (e' 

G (A„(i3^))' /* ^{s, [(1 + (AX,)2 + (An)2) \\X, - Y,\\l 1 



G 



7W 



YWll 



n\nr. 



ds 



ds, 
(4.2) 
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where 

If w G rjjv then, by Lemma l473t we have that 

1 + (AX,)2 + (Ay,)2 < Cn. 

Set 



(4.3) 



Multiplying Equation (14.21) by (p{s,t) and integrating, yields 
VN{t) < Cn [(A„(S^))' + 11/ ip{s,t)VN{s)ds 



(4.4) 











drds. 



By the bounded convergence theorem, we have that almost surely 

ft 



VN{t) 



ds < oo 



and 



drds 0, 



as iV tends to infinity. Then, there exists a random variable N* e N such that 
C cpis, t) ifiir, s)E^ [\\Xr - YrWl hi\n.] drds < iv]v(i), (4.5) 
for aU N > N*. Substituting (g^]) into (031) yields 

r 2 1 f* 

VN{t)<CN (A„(S^)) +1 / ^{s,t)VN{s)ds, 

for all N > N*. Applying now the Gronwall-type Lemma 7.6 in [TS], we deduce 
that VN{t) — for aU N > N* almost surely. Hence, 

P[Xt^Yt, yte [0,T]] = 1, 

and the pathwisc uniqueness property holds. ■ 



5 Existence of solutions 

Let us now introduce the Euler approximations for Equation (|1.2p . Consider 
the framework {Q, P), {J^t, t e [0, T]}, (Xq, B" , W) introduced in Section 2. 
Fix a sequence of partitions 

= < i!^ < . . . < < . . . < = T 

of [0, T] such that 

sup \t^,-t^\^0 

0<j<n-l 
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as n 00. Define X'^{t) = Xq and for n > 1, 



X"(i)=Xo+ / 6(A:„(s),X"(fc„(5)))ds+ / awiknis),X'\k„is)))dWs 



+ / aH(fc„(s),X"(fc„(s)))dBf , 
Jo 



(5.1) 



where 



if t G [i", We will show the following result. 

Proposition 5.1 For any integer N > 1 there exists a random variable Rn > 
0, depending on Xq and , such that, almost surely, 



\xi' - x: 



N 



<Rn \t~s\" , 



(5.2) 



for all s,t e [0, T] and n eN. 



Proof. The proof will be done in two steps. 

Step 1.- We begin by proving that there is a random variable Kp] > such 



that 



E 



w 



\X. 



< Kn, 



(5.3) 



for all t G [0, T] and for all iV G N. 

Note that the paths of X"(fc„(-)) are piecewise constant and the integrals in 
(|5.ip are just finite sums. In the following computations, Cat denotes a positive 
constant that depends on N and the other parameters of the problem, and may 
vary from line to line. From (|5.ip . we have that 



E 



w 



IX 



*llf 



<CNl\Xor+E^ 



t 2JV 

6(fc„(s),X"(fc„(s)))ds 



E 



w 



-E 



w 



cHv(A:„(s),X"(fc„(s)))dW^, 



^H(fc„(s),X"(A:„(s)))dBf 





2N- 



2N- 



2N 



Ai + A2 + Ai 



Using the estimate p.Sp and Holder's inequality, we obtain 



Ai < CnE 
< CnE^ 



\XHK,{s))\ 
(t - sY' 



ds + l 



2N- 



\X^{k„is))fds 



|X"(fc„(s))|'^ds 



Cn 



Cn- 
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We have also that 

A2 < Cn^^ 







■Cat] 



7W 



t 2N- 



2N- 



\a+l 



{t-sf 
= All +yli2. 

Applying Burkholder and Holder inequalities, we have that 



ds 



All < Cat] 



|aw(fcnGs),X"(fc„(,s)))|2^ds 



< Cat / 1 + E 



\X'\kJs))\ 



2N 



ds, 



where we have used the linear growth assumption in (Hctvk)- For the second 
term we have, by Holder and Burkholder inequalities, that 



A12 < Cn] 



1 



2N-1 



{t~sy 



'-{-+ 




1 1/2 + c N 

2 2N ) 

2N 



ds 



Jl<Jw{kn{r),X^{k,,{r))dWr 



it - s)^+4+" 



-ds 



< Cn I {t-s)-^-"E^ 



|fTv^(fc„(r),X"(fc„(r)))|'^dr 



ds. 



Applying now Fubini's theorem and using the growth assumption in (Haw), we 
obtain 



Ai2<Cn / (t-r)-^-"(l + E^ [|X"(fc„(r))|'^]) 



dr. 



Therefore, 



Applying p.7p . we have that 



A2<Cn I (i-s)-^-"E^ [|X"(fc„(s))|'^' 



ds + a 



TV- 



A3 < Ca.A„(B^)2^ 



ds 



By Holder's inequality and the assumptions in (Hcr^f ), we have 



A3 < C'NA^iB^f^ / {{t - s)-2" + 



1 +E 



w 



||A"(fc„(s) 



l\2N 



dr. 
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Putting together all the estimates obtained for A\^A2 and ^3, we obtain 



E 



w 



|2W 



<CA.|Xor^+C^[A„(S^r + 1] 



X \{t-s 



4-" , ^-o^^w 



2N 



ds. 



(5.4) 



Therefore, since the right-hand side of Equation (j5.4p is an increasing function 
of i, we have 



sup E 

0<s<t 



w 



„2N 



< CN\Xor +Cn [A„(B 



H\2N 



t ~ S 



sup it 

0<u<s 



2N 



-1] 

ds. 



As a consequence, by the Gronwall-type lemma (Lemma 7.6 in jl9j). we deduce 
the desired estimate. 

Step 2.- Now we show that there is a random variable Rn such that (|5.2p 
holds. In the sequel, Rn denotes a positive random variable. We have 

2N-\ 



E 



w 



<Cn<E 



w 



6(fc„(u),X"(fc„(u)))du 



+ E 



w 



t 2Ar- 

aM/(fc„(u),X"(A:„(u)))dVK„ 



-E 



w 



aH{K{u),X^{kn{u)))dB^ 



2N- 



= Bi + B2 + B3. 

Applying Holder's inequality, the growth assumption in (H6) and ()5.3p . we have 
that 



\b{kn{u),X'\kn{u))\ 



2N 



Bi < CN{t-sf^-^ 
< i?jv(t-s)2^. 

By the Holder and Burkholder inequalities and using (15. 3|) . we obtain 



B2 < CN{t-s)^-^E^ 



\cTwikn{u),X^{kn{u))\^^du 



<RNit~s) 



N 



Finally, using the estimate (j3.4p and the Holder inequality, we have 



t 2N 

H\2N u „-i27V(1-q)+2q-1 



<Cn^c.{B''Y'' {t-sY 



(r-s) 

Applying this estimate, the assumptions (Rcjh) and ()5.3|) . we obtain 



dr. 



Bz <E'^ K{B"Y'^ (t-s) 

<Rj, (t-s)2^(^-")+^"-^E^ 
<RNit-sf, 



^2N{l~a)+2a-l 



*||aH(fc„(r),X"(fc„(r)))||f 



l + ||A:"(fc„(r) 



(r - sY 



{r-s) 

:2N 



2a 



dr 
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which concludes the proof. ■ 

As a consequence of Proposition l5.ll we estabhsh the tightness of the law of 
the sequence {^"}„gpj in the space Cg of 77- Holder continuous functions, with 
r] < ^, such that 

\f{t)~f{s)\ . 
lim sup — — 0. 

'^-O 0<|t-s|<e [t-sP 

These spaces are complete and separable \10\ . 

Proposition 5.2 Let P" = P o X^, n > 0, be the sequence of probability 
measures induced by X" on Cq . Then this sequence is tight. 

Proof. Fix e > and rj < ^. It suffices to show that there exists a compact 
set K in Cq such that sup„>Q P{X" G K'^) < e. Choose an integer N such that 



2N 



> rj. Let M > be such that 



P[Rn > M) < -. 



(5.5) 



Define a new probability by 

Q{B) 



P{Bn{RN < M}) 



P{Rn < M) 



Then, Proposition 15.11 implies that 



,2N 



{i?,N<M} 



P{Rn < M) 



< MP{Rn < M)-^ \t-s 



N 



By the tightness criterion established in |14j, the sequence of probabilities 
Q o 71 > 0, is tight in Cq . Therefore, there exists a compact subset K in 

Cq such that 

supg(X" e if") < P(i?Ar < Af)~iJ. (5.6) 

n>0 2 

Finally, from (|5.5p and (|5.6p we obtain 

P(X" e K"") < P(X" G Rn < M) + P{Rn > M) < s, 

which allows us to conclude the proof. ■ 

Now we can show the existence of a weak solution for Equation (|1.2p . 

Theorem 5.3 Assume that the coefficients 6, crw and an satisfy the assump- 
tions (Hb), (Haw) o,nd (Han)- If ^ ^ H <a < min|^,/3,|}, then there exists 
a unique weak solution X of Equation 

Proof. The proof will be done in several steps. 

Step 1.- By the Prohorov theorem, the sequence {P", n > 0} is weakly rela- 
tively compact in Cg and exists a subsequence, that we still denote by P", which 
is weakly convergent to some probability P°°. By the Skorokhod representa- 
tion theorem, there exists a sequence of processes {y".P",M^",0 <n< 00}, 
defined on some probability space {il.,T,P) and with values in Cg, such that, 
for every < n < 00, the process (1^", P", W^) has law P" and 



lim ||y"-r= 



IP" - p= 
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almost surely. 

Since, for every n, the process {Y", S", VT") has the same law as (X", , W), 
if we introduce the filtrations 

^<7{Y°°{s),B^{s),W°°{s), s<t}, 

the process W" (resp. W°°) is an (resp. r-dimensional standard 

Brownian motion. Moreover, i?" and B°° are fractional Brownian motions. 

Step 2.- By an adaptation of a result in [23] (page 32) or Lemma 3.1 in [9], 
for any continuous function f{t, x) which satisfies the linear growth property in 
the variable we have that 



lim 

n >oo 



/(fc„(s),y"(fc„(,s)))d,s= / /(,s,y-(s))ds. 



lim 

n >OG 



f{K{s),Y^ (fc„(s)))dw^; 



f[s,Y^[s))dW, 



oo 
s •> 



as n tends to infinity, in probability, and uniformly in i G [0,T]. We have 
also a similar result for the convergence of integrals with respect to fractional 
Brownian motions: 



lim 

n ^oo 



C aH{kn{s).Y"{kn{s)))dB: ^ f cjH{s.Y^{s))dBT, (5.7) 
Jo 



as n tends to infinity, uniformly in i G [0, T] and P-a.s. Let us show the 
convergence (|5.7p . By the linearity of the generalized Stieltjes integral, it is 
clear that 



aH{K{s),Y'\K{s)))dB: 



f aH{s.Y^\s))dBl 
Jo 



< A1+A2, 



where 



and 



<jHikn (5),y"(fc„(s)))d(s^'-sr) 



A. = 



[aniK, {s) , y"(fc„ (s))) - anis, F°°(s))] dB^^ 
Using the estimate p.Sp . we have that 

/* [anis, Y^is)) - (s, Y°°{s))] dB^ < 

Jo a 
< Ck^{B^) 

{{t - sr"' + s-^) [(1 + Ay"(5) + Ar°°(5)) \\Y^\s) - Y°-{s)U ds 



< CK^{B^) ||y" - Y^i^ 1 + ||y"||; + ||y 
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as n tends to infinity,P-a.s. Using the estimate 
(HcT/f ), we have 



and the assumptions in 



< CA„(i?°°) \\aH{kn (s) , r"(A:„ (s))) - gh{s, Y'\s))\\^^^ 

< CA„(i3-) Wanikn {■) ,Y^\k,, (•))) - cTHi;Y^\-))\\lo ih +h + h) , 
where e > is a small positive number that depends on a and /? and we have 

h = hnikn (•) , y"(A;„ (•))) - cTff (•, r"(-))|lL"^ > 



h = 



T 
"'0 
<Ci+C2 



aH(s,y"(s))-aH(r,y"W)l 



(s-r) 

|y" (s) - y"(r) 



-drds 



<Ci + C2 iiy 



n||(l-e)'? 



(s-r)' 



-drds 



where we have used the Holder continuity in time in assumption (Kan), and 



7^ 

"'0 



anikn (s) , y" (fc„ (s))) - t7H(fc„ (r) , y"(fc„ (r))) 



< c 



T ps 



^0 



(s-r) 

|fcn (s) - fcn + |A:„ (s) - fc„ (r)!'!-^)" 



-drds 



(s-r) 



Q + l 



drds 



We can compute this last integral, using the partition on the interval and de- 
composing the integrals in finite sums. This integral is uniformly bounded in n. 
Therefore 



(aniK, (s) , y"(fc„ (s))) - anis, y"(s)))di?= 



0, 



as n tends to infinity, P-a.s. In order to show the convergence of the term Ai, 
we use again the estimate p.4p and Lemmas 7.4 and 7.5 in [19]. We obtain that 



< ||aH(fc„(s),y"(fc„(s)))iu^,A„ (sr-sD 

<ciii?"-i3-ir^, 

where e > is a small positive constant which depends on a. Therefore, Ai 
converges to zero as n tends to infinity, P-a.s. 

Step 3.- Recall from step 1 that (y",B",W^") and (X",P^,iy) have the 
same laws. Moreover, W" is a standard Brownian motion in the appropriate 
filtration and P" is a fractional Brownian motion. Therefore, our processes 
satisfy the stochastic differential equations 



V = ^o" + / b(kn{s),Y'\k,,{s)))As 



+ / aw{kn{s),Y^{kn{s)))AW: 



aH{k^{s),Y^{k^{s)))dB: 
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almost surely. So, by step 2, when n tends to infinity, we obtain 



Jo Jo Jo 

Therefore, Y°° satisfies p.2p with the driving noises W°° and B°° . 

The sample paths of Y°° belong to Cq C Wq'°° almost surely, and further- 
more, by (|5.3p . we have that 



T 







ds < oo. 



Therefore, by Definition [HI {Y^,W°^,B°°) is a weak solution of (H^. 
We can now proceed with the proof of Theorem [ 



Proof of Theorem 12.21 The uniqueness is a consequence of the general 
pathwise uniqueness proved in Theorem 14.41 For the existence of a strong so- 
lution we can make use of the classical result by Yamada and Watanabe [21], 
which asserts that pathwise uniqueness and existence of weak solutions imply 
the existence of a strong solution. The main difference with the classical proof 
is that here we have two random sources independent of the Wiener process W, 
the initial condition Xo and the fractional Brownian motion . It suffices to 
replace R"^ by the product space R'^ xC([0, T])™, endowed with the product 
measure jiy. v, where fj, is the law of Xq and v is the law of B^ on the space of 
continuous functions. ■ 
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